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ABSTRACT 

The step height effect of radially inward expanded air flow stream on the heat transfer process in a concentric 
circular annular passage is studied experimentally. The airflow separation, reattachment and redevelopment took place in 
the test section which its inner tube subjected to a constant heat flux boundary condition. The experimental apparatus 
comprised of concentric tubes to form an annular passage with a sudden enlargement in the passage cross-section created 
by increasing the outer tube diameter of the annular inner tube at the entrance section prior to the test section. The inner 
tube of the test section was made of aluminium having 25 mm outside diameter and 350 mm heated lengths, which was 
subjected to a constant wall heat flux boundary condition.The investigation was performed to cover a Re range of 3000 - 
11000, inner tube heat flux varied from 800 W/m 2 to 1750 W/m 2 and the enhancement of step heights were S = 0 ( without 
step), 6.5 mm, 12.5 mm, and 18.5 mm which refer to Dfd = 1, 1.5, 2 and 2.5, respectively. 

Results reveal for all cases, that the local heat transfer coefficient increased against increasing the heat flux 
and/or Re. The step height has an eminent effect on heat transfer in the separation region which decreases with the raise of 
the step height. The passage without step (S=0) results were correlated and compared with the availableturbulentforced 
convectionentry length correlations for same geometry and boundary conditions. The step height effects show, for S= 6.5 
mm a 25% on the average improvement in heat transfer for a full Reynolds number tested while for another two step 
heights the results show a limited improvement in low Reynolds number to a 16.5% reduction in heat transfer due to 
separation zone effect. The velocity profiles demonstrate that the position of reattached point behind the step increases with 
the rise of step height. The present results show a good agreement with the available previous works and have followed 
similar trends. 
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INTRODUCTION 

The separation of boundary layer takes place when the part of the boundary layer close to the wall reverse in flow 
direction and that occurs as the overall boundary layer thickens suddenly to force it off the wall by the reverse flow. For 
internal flow, such as a rapid expansion of pipe, an adverse pressure gradient met as the flow expands creating an extended 
region of separated flow. Separation generates two parts; the recirculating flow and dividing streamline at the central 
region of the pipe. The point at which the dividing streamline attached the wall is called the reattachment point. As the flow 
goes further downstream, it eventually redevelops and settles down to achieve an equilibrium state [1]. The flow separation 
as one of the viscous flow problems is very important for science and industry. As per classical concept, flow separation is 
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caused due to viscosity. Therefore, it is often expressed as “boundary layer separation”. Generally there are two types of 
flow separation, separation of external and internal flow [2]. The classification of separation can also sort related to flow 
behaviour as separation with or without reattachment. Typical examples of the former are the flow over a backward facing 
step, the leading - edge separation bubbles of air foils and blunt plates or blunt circular cylinders. The separation flow 
described as reattaching when the interaction exists between the vortices and the solid surface and it described as without 
reattaching when the interaction is between vortices shed from the separation point [3]. In general, this phenomenon is 
encountered in some engineering applications, such as flow over air foils at large angles of attack, in channels where area 
suddenly increases, in wide angle diffusers [4-9], heat exchangers, combustors, nuclear reactor cooling channels [10], in 
power plants [11], gas turbines and electronic circuit panels [12-13]. 

Forced convection heat transfer in a concentric circular annular passage has been studied experimentally and 
analytically by Kays and Leung [14]. The experimental solutions established for air flow with thermally entry length and 
second kind boundary conditions and employed in a broad range of Reynolds numbers with either outer tube heated or 
inner tube heated with constant heat flux. The experimental data converted to an empirical equation which was in a good 
agreement with analytical solutions. A correlation equation has been presented for air flow in annular passage with constant 
heat flux subjected to the inner tube only and for annulus radius ratio r* = 0.2 (r*=r/r 0 ) for as asymptotic Nusselt number 
as: 

Nu Dh = 0. 0324 Re^ 62 '. (1) 

Using THMOD2 code, Bsebsu and Bede [15] presented a theoretical investigation dealing with a down flow 
single-phase forced convection in a narrow vertical annuli sub channel (WWR-M2 channel). The investigation main 
objective was to examine the applicability of available turbulent heat transfer equations in a narrow vertical annuli channel. 
The module represents a 3 mm spacing (gap) and 600 mm active length in fuel elements used in thermal hydraulic analysis 
tasks of the WWR-M2 research reactor. The investigation concluded that by using the equivalent hydraulic diameter the 
existing correlations are applicable to WWR-M2 channel. A new equation has also been developed for sub channels of 
WWR-M2 for either one or both side heated conditions. 

The heat transfer characteristic in a double pipe with 560 pin fin heat exchanger was studied theoretically by 
Nieckele and Saboya[16]in which the air is flowing in the annular channel while the water is flowing in the inner circular 
pipe. The results depicted in dimensionless form and in terms of average Nusselt number and friction factors as a function 
of the Reynolds number ranging from 13,000 to 80,000. A comparison with the available smooth sections was presented. 
Filetti and Kays [17] experimentally studied the local heat transfer rates to air flow close to the entrance of a flat duct in 
which there is an abrupt symmetrical enlargement at the flow cross section. Reynolds numbers varied from 70,000 to 
205,000 and two enlargement area ratios (2.125 and 3.1) were considered.lt was found that such a flow was characterized 
by a long stall on one side and a short stall on the other. Maximum heat transfer occurs in both cases at the reattachment, 
followed by decay towards the values for fully developed flow. The results also revealed that the Nusselt number along the 
annulus inner tube downstream the step for the separation, reattachment and redevelopment regions has the same behavior 
as described by Vogel and Eaton [12]. Heat transfer in the separation region behind an obstacle was studied experimentally 
by Davletshin et al. [18]fora wide frequency range of superimposed free stream pulsation. The heat transfer coefficient was 
evaluated by solving an inverse non-stationary heat conduction problem based on experimentally measured, the wall 
transient temperature. The pulsation flow shows substantial heat transfer intensification in the separation region. 
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The effect of Reynolds number and the duct divergence angle on the flow separation characteristic in a conical 
duct was studied numerically by Sparraw et al [19]. The flow considered as fully developed at the diffuser inlet while the 
diffuser exit connected to a long circular pipe. A moderate-turbulence and transition Reynolds number at the diffuser inlet 
accommodated to get a possible laminarization of flows using a universal flow regime model. The investigation outcome 
found that the separation angle for Reynolds number less than 2000 occurs with a diffuser expansion angle of 5 degrees. 
The results conceal a prior guideline that separation first occurs at a divergence angle of 7 degrees. The simulation shows 
that for divergence angle from 20 to 30 degrees, the separation occurs in all investigated Reynolds numbers. Togun et al. 
[20] investigated experimentally the effect of step height on the outward separation air flow in a concentric annulus in 
which the annulus outer surface is subjected to a constant heat flux. Investigation covered in a Re range of 17050-44545, 
heat flux varied from 719 W/m 2 to 2098 W/m 2 and the enhancement of step heights were, S= 0 (without step), 6 mm, 14.5 
mm and 18.5 mm, which refer to d/D = 1, 1.16,1.53 and 1.80, respectively. The effect step height of flow separation has 
been correlated and the results reveal that the maximum enhancement of 18% occurs for step height 18.5 mm. Local heat 
transfer coefficient demonstrates a gradually improving in heat transfer directly after the step in the separation region to 
reach maximum value in the reattached point then it decreases gradually as the boundary layer builds-up in the redeveloped 
region. 

The entry lengths for turbulent flow are typically short in comparison with the laminar flow entry length often just 
10 tube diameters long, and therefore the Nusselt number determined for fully developed turbulent flow can be used 
approximately for the entire tube. This simple approach gives reasonable results for heat transfer for long tubes and 
conservative results for short ones. The correlation for the average heat transfer coefficient in a short annulus is very 
limited and the condition of applying notorious equation of Dittus and Boelter [21] is far from the present experiment 
condition. Therefore, to compare the present experimental result, the empirical equation should be taking into account the 
effect of Reynolds and Prandtl numbers, the position of the passage entrance, and the annulus side heated. Petukhov [22] 
and Gnieliski [23] correlate for circular tube and the available data somewhat better than any other over the range of 
Prandtl number from 0.5 to 2000, and Reynolds number from 2300 to 5xl0 6 . Kim et al. [24] recommended the following 
empirical equation for the low turbulent Reynolds number forced convection heat transfer in the entrance of the annulus: 

Nu = 0. 039 C Therm Re° 74 Pr 04 . (2) 


WhereC X | lerm is a function of Re and the normalized distance from the start of the heating zone, x/D h , where D h 
represents the hydraulic diameter: 


C Thm = 1.0 + exp (o. 08. (£))(£)' 


0.65 + 


( 2925 
V Re 


)© _ “ 


(3) 


Low Reynolds turbulence model results are compared with the experimental data available and the predicted by 
v 2 -f model and the comparison was found much closed. A correction to the above forced convection equation is suggested 
by Petukhov and Roisen [25] for the annulus with inner tube wall heated with constant heat flux while the outer wall is 
insulated: 


^ = 0.86 
Nu F 




(4) 


Where Nu P represents the output of the equation presented by Petukhov [22] or Gnieliski [23]for circular tube 
while Nu c is the correctedNusselt number value for an annulus has the outside diameter D 0 and the inside diameter Dj. 
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There is a shortfall of experimental observation to describe the separation and heat transfer processes in annular 
passage especially for short passage. The motivation of present investigation is to add a modest contribution to the flow 
separation phenomenon. This work focuses to study the effect of flow separation of the heat transfer process for an 
axisymmetric, transient to turbulent separated flow in an annular passage. The air flow separation was induced with 
different heights of circular steps by increasing the inner annular tube diameter at upstream of the heated test section. 
Reattachment and redevelopment of the flow occurred behind the point of separation. The test pipe was heated electrically 
to conduct the experiments at constant heat flux and the degree of separation is varied by using different tube diameters 
while the Reynolds number was maintained in the same range by regulating the gate at the entry of the blower. 

Experimental Apparatus 

The apparatus was designed and constructed to have an axisymmetric separated and reattached air flow in an 
annular passage. The full set-up with annular shaped test section is shown in Figure 1. The passage proceeds upstream with 
a variable entrance calming section. The apparatus can provide different air velocities (variable Reynolds numbers) and 
heat fluxes at the test section. In general, the apparatus consists of two parts, the open air circuit and the measuring system 
(Table 1). 

An open-air circuit, as shown in Figure 2, consists of a small centrifugal fan (A), delivery section (B), settling 
chamber (C) and the circular concentric annular tubular passage (D, E and F). The annular passage has started with an 
entrance section (D) and has ended with an exit section (F). The air delivered by a variable speed centrifugal fan (A) and 
was controlled by a fan speed regulator (Variac) and fine adjusting the position of the lever of the inlet gate valve (G) at the 
fan entrance. The velocity of the air was measured by using a hot wire anemometer (H) located at the end of thecircular 
delivery section (B). Air then passes through the settling and the calming section (C) which equips with two layers of fine 
gauze (I) and a layer of the flow straightener (J) to streamline the incoming air before delivering it to the passage bell 
mouth (K). Then, the air passes through the circular concentric tube annular passage (D to F). It is constructed of an 
unheated aluminum outer tube (Q) of constant inner diameter of (140 mm) and elongated to cover all the annular passage 
while the passage inner part has two tubes with variable diameter tube (M) located at annulus calming entrance section and 
the heated test tube (L) with constant diameter of (25 mm) located at annulus passage test section. Both tubes held and 
supported horizontally and centrally with the help of the threaded rod (N) which rested in the exit centered Teflon piece 
(R2) and fasten withhexagonal nut (02) from one side while the other side rested on the cross metal piece (P) in the 
settling section and fasten with hexagonal nut (Ol). This combination allows the passage inner tube to be assembled and 
dismantled easily for changing the entrance section inner tube and to adjust its position horizontally and centrally in the 
passage. The calming entrance section inner tube outer diameter varies according to the separation step height case study. It 
has four variable diameters D, of 25, 37.5, 50 and 62.5 mm which provides a separation diameter ratio (D/d) of 1.0 (with 
no separation), 1.5, 2.0 and 2.5, respectively. The combination of annulus inner tube (M) and the two Teflon pieces (Rl) 
(connecting the two tubes (M) and (L) and on which separation takes place) and (R3) (adopted with the bell-mouth (K) to 
minimize the flow disturbance at the annular passage entrance) is sliding on the threaded rod during the assembly. The 
annulus outer tube (Q) in the test section part has seven holes (W) which allow the Pitot static tube to be penetrated during 
the velocity profile test. The Pitot Static tube is connected with the micro-manometer via a flexible connection (U). To 
measure the bulk air temperature at the passage inlet and outlet, four thermocouples have been inserted in the flow passage 
one thermocouple (SI) to measure the air inlet temperature and three thermocouples (S2) to measure the bulk air 
temperature at the passage exit. 
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The air at the test section E, (which is the last part of the open air circuit) experienced a separation, reattachment 
and redevelopment at the passage downstream. It was heated by the annulus inner tube of that test section (L) at constant 
heat flux (CHF). The outer surface of connecting Teflon pieces (R2) was especially machined to be flush mounted with 
both calming section inner tubes (M) and test section inner tube (L) to prevent any flow disturbance at the separation point. 
Teflon was chosen for easy machinability and low thermal conductivity in order to reduce the test section end losses. All 
pieces were machined from a Teflon rod (75 mm) in diameter. Two layers of fibre glass (40 mm thickness) were used to 
cover the outside surface of the annular passage outer tube (Q) to reduce the heat losses radially outwards the annular 
passage. The fan discharges the air to the system via a flexible connection (V) to prevent fan vibration to be transmitted to 
the system. 

Heating System 

The test tube was heated electrically by using an electrical heater as shown in Figure 3. The heater consists of a 
nickel-chrome (Nickhrom) wire (A) uniformly wounded to produce a constant heat rate per unit length. The heater was 
isolated electrically from the aluminium tube (B) by winding with a very thin asbestos tape and the wire itself was isolated 
by fine ceramic beads. A 10 mm wide tape spacer has been wounded while winding the heater to assure uniform heater 
winding. The temperature along the outer tube surface was measured by 16 thermocouples (D) located along the tube, 
starting from the separation point (X = 0) up to the end of the test section as shown in Figure 3. At first, thermocouple 
positions were marked along the entire tube surface, and then 1.5 mm holes were drilled on the wall at 2 mm depth. The 
drilled holes were then chamfered by a 2 mm drill and cleaned properly. Thermocouples hot junctions (E) were inserted 
inside the tube and installed and secured in cleaned holes with adhesive (F). The outer surface of the test tube then cleaned 
to remove the excess adhesive and to make the measuring point flash mounted with the tube surface.To avoid convection 
currents in the space between the heater and the cylinder inner surface, this space is filled and compacted with magnesium 
oxides (C). Magnesium oxide was selected as it has relatively high thermal and poor electrical conductance. 

To determine the heat loss by conduction from the test tube ends to the Teflon pieces, two thermocouples (G1 and 
G2) inserted on each side of the test tube in the end Teflon pieces (HI and H2) as shown in Figure 3. Knowing the thermal 
conductivity of Teflon and the distance between the thermocouples, the end conduction could be calculated. The maximum 
value of this, varied in the experiments between 1.23% and 2.5% of the total heat input. To evaluate the conduction heat 
loss from the whole fiber glass insulation cover the whole passage, the average temperature of the four thermocouples (I) 
under fiberglass of the test section and the average temperature of two other thermocouples located (2 mm) apart from the 
first two were recorded. Then, by knowing the fiber glass thermal conductivity, the heat losses could be calculated. The 
amount of heat loss varied in the experiments and the maximum value of this was less than 0.25% of the total heat input. 
All thermocouples used in the test rig are K type, 0.1 mm diameter wire and with two-layers of asbestos insulation. All the 
thermocouples were calibrated in the laboratory using the boiling and melting point of pure chemical substances. The 
thermocouples were calibrated without lead and with lead (similar to experiment condition) to determine the correction 
factor for applying leads and the results show no appreciated effect on the calibration. The thermocouples were made by 
fusing the ends of two wires by an electric spark in an oxygen free environment. 

Data Analysis Method 

The following simplified steps were used to evaluate the heat transfer to the flowing air in the expanded annular 
passage while the inner surface of outer tube was subjected to uniform wall heat flux boundary conditions. The total heat 
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transfer can be calculated using: 

Qt = VI . (5) 

The convection and radiation heat transferred from the test section inner tube surface is calculated by: 

Qc+r — Qt — Qcond. ••• (6) 

Where Q con d is the conduction heat losses through Teflon end pieces and Q c+r is the convection and radiation heat 
transfer. The conduction heat transfer from fibre glass can be represented by: 


Qcond ^Teflon^Teflon 


ATreflon 

■ lX Tel hill 


.(7) 


Where AT 


Teflon 


is the temperature difference inside each of the test tube Teflon end pieces, Ax^y^ represents the 


distance separate the thermocouples, k Te f lon is the Teflon thermal conductivity and A Te ^ lon represents the Teflon end piece 
cross-sectional area. 


The convection and radiation heat flux can be represented by: 

q C+ r = ^ . 


( 8 ) 


Where A s — n d. L, L and d represent the test tube length and diameter, respectively. The radiation heat flux can 
be represented by the following Eq.: 

q r = crF 12 (Tt 1 -T 4 a2 ) . .(9) 

Where ois the Stefan-Boltzmann constant (5.669 x 10" S W/K 4 m 2 ); while F is the shape factor is equivalent to: 

1 


F 12 = 




....( 10 ) 


Ei and s 2 are the annulus surface emissivity values (outer surface of the annulus inner tube and the inner surface of 
the annulus outer tube in the test section), T) is the average outer surface of the test section inner tube temperature and T 2 is 
the average inner surface of the test section outer tube temperature. The convection heat flux can be represented by: 

. .( 11 ) 


Re Rc+r Rr 

The convection heat flux is used to calculate the local and average heat transfer coefficients as follows: 
q c 


h x = 


( T SX - T BX) 

All the properties were evaluated at the mean film temperature which evaluated by the following Eq.: 
t sx- t bx 


( 12 ) 


Tfx = z 

Where T BX is the local bulk air temperature and Tfx is the local mean film temperature. 
The local Nusselt number based on the diameter (Nu x ,Dh) can be determined by: 

h x .Dh 


(13) 


Nu 


X,Dh — 


Kf 


(14) 


The local Nusselt number based on the axial distance (Nu x ) can be evaluated by the following Eq.: 
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Nuv 


hxX 

Kf 


■(15) 


The average values of the Nusselt number (Nu L ) can be calculated based on the calculated average heat transfer 
coefficient by: 


h L = 7/0 hx- dx 


(16) 


The average Nusselt number based on the calculated average surface temperature and average bulk air 
temperature are evaluated: 


q c 


OsL - T BL ) 


(17) 


Where 


T L = i/ Q L T x .dX. 

rr _ T X= 1+ T x= 0 

*BL — 2 . 

All the average properties were evaluated at the mean film temperature as: 
T _ (t l ~Tbl) 


The local Reynolds number based on the hydraulic diameter (Re Dh ) can be calculated by: 


Re 


_ Pf 

Dh — — 


.U.Dh 

Pf 


( 18 ) 

(19) 


( 20 ) 


( 21 ) 


And the local Reynolds number based on the distance (Re x ) can be determined by: 


Re 


X 


P/-U-X 

Pf 


( 22 ) 


All the air physical properties p f , p f , v f and k were evaluated at the average mean film temperature (T f ) Holman 
(2010) [8]. A computer program in FORTRAN was written with algorithm capable to manipulate the data set calculation of 
air all properties locally and doing and select best equations which fit surface temperatures and local heat transfer 
coefficients and then doing the numerical integration to obtain the average value for heat transfer coefficient and 
dimensionless groups for the whole cylinder. 


RESULTS AND DISCUSSIONS 


To cover all experimental parameters, check the measurement consistency and to test the velocity profile variation 
along the passage, a total of 70 test runs were conducted. The parameters to cover include a straight annular passage 
without step (S = 0) and three passages with step heights corresponding to S = 6.5mm, 12.5 mm and 18.5 mm hold the 
passage expansion ratios (D 2 -d 2 ) / (D 2 -D , 2 ) = 1, 1.045,1.11 and 1.2 respectively. The test section inner tube heat flux was 
varied from 800W/m 2 to 1750 W/m 2 while Reynolds numbers Re D wasvaried from 3000 to 11000. The velocity profile test 
runs cover the annular passage with steps and for seven axial positions. The data reproducibility is good and remains 
within the confidence level of 95% and average rms error less than 1.5%. 
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Variation of Surface Temperature 

The variation of surface temperature along the test section inner tube which may affected by many variables, such 
as flow velocity, surface heat flux and upstream flow separation created by the step. The temperature variations for a 
selected runs are demonstrated in Figure 4,5 and 6. 

Figure 4 shows the surface temperature variations along the test section heated inner tube for different heat fluxes 
with a specific Re D = 3000 and step height S=18.5 mm i.e. (D/d = 2.5). The general distribution shows a reduction of 
surface temperature immediately behind the step at the test section entrance (x = 0) to reach the minimum value at specific 
axial position at which the flow reattachment expected to occur after the separation region due to the turbulence which 
increased the heat transfer. The surface temperature then increased gradually to reach a maximum value then it decreases at 
the test section exit due to the tube end losses. 

Figure 5 depicts the variation of Re number effect on the temperature distribution of the test section inner tube at 
specific heat flux equal to 1750W/m 2 and step height equal to 18.5 mm. The figure reveals that the surface temperature 
decreases as the Re D increases for the same heat flux and step height. 

Figure 6 presents the temperature distributions for selected runs using the step height as a variable parameter. It is 
detected that while the surface temperature for the straight passage increased gradually for whole the heated surface the 
other three cases show that the surface temperature increases as the step height increases especially in the separation 
region, directly behind the step for the same Reynolds number and surface heat flux. The figure also reveals that the 
position of the minimum surface temperature moves downstream as the step height increases to indicate that the separation 
region behind the step increases as the step height increases. The surface temperature distribution over the separation 
region has a similar trend of a straight annulus case with a gradual increase of surface temperature for the rest of the tube. 
The minimum temperature represents the expected flow reattached point. Separation produces a recalculated flow moving 
towards the step at the heated tube surface which then drifts away from the step by the main separated flow. 

Evaluation of Local Heat Transfer Coefficient (h x ) 

The variation of the local heat transfer coefficient h x along the test section inner tube (X) for selected runs is 
shown in Figure 7 and 8. Figure 7 depicts the effects of Reynolds number variation on the h x distribution for a 
corresponding heat flux equal to 1750 W/m 2 . Thedistribution of local heat transfer coefficient h x increases until it reaches a 
maximum value at the reattachment point (maximum local heat transfer point) after which h x decreases to grow 
asymptotically while moving downstream along the X-axis towards the test section exit. 

Figure 8 shows the variation of local heat transfer coefficient h x with straight annulus (without step) and three 
different expanded passages due to different step heights. The straight annulus shows completely different distribution 
compared with the other cases (with step) as the distribution decreases sharply at the test section entrance then decreases 
gradually to end asymptotically at the section exit. The high h x at x=0 is due to zero thermal boundary layer at the test 
section entrance. For separated flow ( with step) and same Reynolds number and heat flux, the local heat transfer 
coefficient h x increases as the step height increases, especially in the separation region due to still recirculating induced 
vortex which move towards the step at the inner tube surface and drifting away by the separated flow. The still 
recirculating vortex reduces the heat transfer rate directly behind the step. The figure also reveals that h x values for D /d 
equal to 1.5 shows a quite improvement in heat transfer rate in comparison with straight annular passage and that attributed 
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to a limited and unclear separation zone for this step height and this step height also increases the case turbulent intensity. 
The other two cases (D, /d equal 2.0and 2.5) show different behaviour as in the separation region, the decline of h x is quite 
clear in comparison to h x obtained for the flow in a straight annular passage. 

Flow Profiles 

The variation of velocity profiles for the separated flow along the test section annular passage and for three step 
heights 6.5 mm, 12.5 mm and 18.5 mm are presented in Figs. 9, 10 and 11, respectively. The three figures are presented for 
the same Reynolds number equal to 5000 and same heat flux equal to 1750 W/m 2 Each figure comprises a seven profiles at 
X equal to 0 mm (at the test section entrance), 25 mm, 50 mm, 75 mm, 125 mm, 175 mm and 225 mm. The profile shows 
the variation of dimensionless axial velocity U against the dimensionless radial position R when r measured from the 
annulus center line. 

Figure 9 reveals that all the profiles in this figure have covered the entire passage with the profile at the test 
section entrance (X/D h = 0) shows that the flow could be regarded as a fully developed turbulent annulus flow for Re > 
2000. The lower part of the profile at X/D h equal 0.215 just behind the step is sharply retreat and the first radial position 
near the tube surface is hardly to measure as the micro-manometer recorder shows unsettled reading with a lot of 
fluctuation due to higher expected turbulence intensity. Due to expected symmetry, the other profiles show a clear 
progressive development in the internal part annular boundary layer adjacent to the test section inner tube. The flow 
behaviour combined with the result of the temperature distribution reveals a very limited separation zone for this step 
height. 

For step height D t /d. Figure 10 also shows the profile at the test section entrance (X/D h = 0) could be regarded as 
a fully developed turbulent annulus flow. The lower part of the first profile at X/Dh equal 0.215 just behind the step is also 
sharply retreat and the first radial position near the inner tube surface is fluctuating around zero and hardly to measure 
which indicates that the flow in this goes in the opposite direction and it is not a part of the main flow.The spatial growth of 
velocity fluctuations from the measurements demonstrated a tendency that casting a vortex behind the step. Therefore, this 
profile is left uncompleted as the lower part of it is dominated by recirculating vortex in the flow separation zone. The 
other profiles show clear progressive development of the boundary layer adjacent to test section inner tube. The 
reattachment length could be extended between axial local position X equal to 25mm and 50 mm. 

Figure 11 shows the velocity profiles variation along the test section for step height of 18.5 mm (D, /d = 2.5) 
which reveals that the profile at the test section entrance (X/D h = 0) could be regarded as a fully developed turbulent 
annulus flow and the reattachment length increases as the step height increases as the second and third profiles suffer a 
sharp retreat while first two local radial positions in the second profile and first local position in the third profile near the 
inner tube surface is fluctuating around zero and showing a reverse flow to the main passage flow. Again the other profiles 
show a progressive development of the annular inner boundary layer. The reattachment length behind the step could be 
extended between the axial local position X equal to 50mm and 75 mm. 

Development of Average Heat Transfer Correlation 

To estimate the effect of step height on the present experimental data, it is required first to correlate the straight 
circular annular passage data (without step) and compares it with the available equations dealing with the present geometry 
and boundary conditions.The correlation for the average heat transfer coefficient in a short annulus and annulus entry 
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length is very limited and the conditions of applying Dittus and Boelter [21] equation (fully developed flow) are far from 
the present experimental condition. All the references indicated that the result of heat transfer at the entrance region 
(developing flow) is higher than the fully developed flow for both laminar and turbulent flows. The result obtained by 
applying Petukhov [22] and Gnieliski [23] (for Re>2500) correlations and the correction for annular passage Eq. (4) 
suggestedby Petukhov and Roisen [25] is approximately equal to the Eq. (2) recommended byKim et al. [24]. Figure 12 
shows the present correlation for the present straight annulus with inner tube subjected to constant heat flux is compared 
with both Kim et al. [24] Eq.(2) and the well-known fully developed flow proposed by Dittus and Boelter [21] equation. 
The figure shows that the results of the present correlation are higher than the results obtained by applying Eq. (2) in the 
range from 5% to 14% for the whole applied Reynolds number. 

The relationships between the average Nusselt number (Nu LDh ) and the average Reynolds number (Re LDh ) are 
shown in Figs. 12-13 to represent the measured heat transfer along the whole test section tube (separation, reattachment 
and redevelopment regions) for different step heights. The whole range of Reynolds numbers and heat fluxes are covered. 
The coefficients of these correlations are presented in Table 2, and the general correlation equation is presented in Eq. (23). 

Nu LDh = C Re" Dh . .(23) 

For step height 6.5mm (Di /d = 1.5), the correlation shows an improvement in heat transfer process for full 
experimental range of Reynolds and it varied from 41.5% to 9.7% as Reynolds number changed from 3000 to 11000. The 
correlations for another two step heights show a limited improvement in heat transfer at low Reynolds which at Reynolds 
number equal to 3000 equivalent to 12.4% and 4.08% at step height 12.5 mm and 18.5 mm, respectively. For the Reynolds 
number 11000 both steps show a reduction in heat transfer which equivalent to -16.84% and -16.34%at step height 12.5 
mm and 18.5 mm, respectively. 

Figure 14 shows a correlation between the local maximum Nusslt numbers based on D h divided by dimensionless 
parameter (D h /S) with the local Reynolds number based on D h . The correlation shows the relation between the maximum 
local heat transfer coefficient h x , maX j step height S and the Reynolds number based on D h . The linear fit of the experimental 
data for the three steps is presented by the following equation: 

Nu maX ' Dh = 0. 0346 Re° D 6 h 19 (^) . .(24) 

CONCLUSIONS 

From the experimental investigation into the range of Reynolds numbers, heat fluxes and step heights, the following 
conclusions can be drawn: 

• In comparison with the straight passage (without step), the surface temperature distribution for step annular flow 
shows a notable increase in magnitude for same Re number and heat flux and the increase rises with the increase 
of step height. 

• Increase the step height and the flow causes the surface temperature coming down up to a specific point (lowest 
temperature) along the test section inner tube and it then increases. The lowest temperature represents the flow 
reattachment point. The position of the lowest temperature is mainly a function of step height. 

• The local heat transfer coefficient (h x ) increases as the Reynolds number increases for all cases with or without 
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step. In the separation region (recirculation zone) the local heat transfer coefficient improves until it reaches a 
maximum value (h Xjmax ) at the flow reattachment point. 

• Reynolds number variation has limited effect on the flow reattach point in comparison with the step height and an 
increase of the flow Reynolds number moves the reattachment point downstream away from the step. The location 
of maximum heat transfer coefficient moves downstream as the step height increases. 

• Flow velocity profiles reveal that the low step heighttested (S = 6.5 mm) does not demonstrate a clear separation 
zone with a clear recirculation flow behind the step in comparison with the other high step height cases (S = 12.5 
mm and S = 18.5 mm) 

• Average Nu number and Re number can be correlated for all cases (with or without step) of flow configuration in 
the annular space. The maximum Nu number, step height and Re number can also be correlated for all with step 
cases. 
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Table 1: Specification of Equipment and Measuring Instruments 


Instrument 

Specifications 

Micro-manometer 

Furness Controls ltd. Type FC014, Model 1, total resolution 10 1:20,000, accuracy 
±.0004"H20, 12” Pitot static probe and connections. 

Hot-wire anemometer 

E+E elecktronic,EE75, Model B,0-40 m/s, +0.03 m/s 

Analog Ammeter: 

Goerz, Type 324760, No. 691907 

Centrifugal Fan 

Marquis Type, : 30 phase, 6/7 (m3/min), power 0.3 kW 

Digital Voltmeter 

Philips, Type PM2404, +0.83% full scale 

Digital Electronic 
Thermometer 

Comark, Type No. 6110, serial No. 11725, Resolution 0.1 C 

Heater 

Nickel Chrome wire (NIC80-040-250) 1 mm Dia. 

Transformer (Variac) 

Lubcke, Type R52, No. 531.51 

Selector Switch 

Comark Type No. 1695, serial No. 24054& 

Comark Type No. 1695, serial No. 11725 

Thermocouples wires: 

Two layers Glass covered Iron-constantan Type k 0.3 mm Dia. 


Table 2: Values of S, C and Slop n at Different Step Height 


Diameter Ratio (Dj/d) 

Step Height S mm 

C Value 

Slop n 

1 

0 

0.197 

0.664 

1.5 

6.5 

1.339 

0.468 

2 

12.5 

1.419 

0.432 

2.5 

18.5 

0.787 

0.496 
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Figure 1: Annular Passage Configuration and Dimensions 
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Figure 2: Schematic Diagram of the Experimental Apparatus 
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Figure 3: Test Section Heating and Measuring Systems Details 



Figure 4: Variation of Local Temperature T x along the Test Section Inner Tube, 
(for Re = 3000 and D/d = 2.5) 
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Figure 5: Variation of Local Temperature T x along the Test Section Inner Tube, 
(for q = 1500 W/m 2 and D/d = 2.5) 
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Figure 6: Variation of Local Temperature along the Test Section Inner Tube, 
(for Re Dh = 11000 and q = 1750 W/m 2 ). 
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Figure 7: Variation of Local Heat Transfer Coefficient along the Annulus Inner Tube 

(for q = 1750 and D/d = 2.5) 
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Figure 8: Variation of Local Heat Transfer Coefficient h x along the Annulus Inner Tube 

(for Re = 11000 and q = 1500 W/m 2 ). 
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Figure 9: Variation of Pitot Static Measurements of Axial Velocity Profiles along the Test Section Annular Passage 

(for Re = 5000, q = 1750 W/m 2 and D/d = 1.5) 



Figure 10: Variation of Pitot Static Measurements of Axial Velocity Profiles along the Test Section Annular Passage 

(for Re = 5000, q = 1750 W/m 2 and D/d = 2) 
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Figure 11: Variation of Pitot Static Measurements of Axial Velocity Profiles along the Test Section Annular Passage 

(for Re = 5000, q = 1750 W/m 2 and D/d = 2.5) 



Figure 12: Comparison of Present Straight Annular Correlating (without Step) with Ref [24] and Dittus and Boelter 

[21] for Fully Developed Flow Equation 
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Figure 13: Comparison of Correlating between Average Nusselt Number with Reynolds Number for All Steps and 

the without Step Annular Passage 
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NOMENCLATURE 


Symbol 

Description 

Units 

A s 

Test section inner tubeoutersurface area 

[m 2 ] 

c P 

Specific heat 

[KJ/Kg.°K] 

d 

Test section inner tube outer diameter 

[ml 

Dh 

Test sectionhydraulic diameter (Dh = (D-d)) 

[ml 

Dh 

Hydraulic diameter at entrance section (D t ,=(D-D i )) 

[ml 

D 

Whole annulus passage outer tube inner diameter 

[ml 

Di 

Entrance section inner tube outer diameter (variable) 

[ml 

F 

Annulus radiation shape factor 


h 

Heat transfer coefficient 

[W/m 2 .°Kl 

I 

Heater current 

[Ampl 

k 

Thermal conductivity 

[W/m.°Kl 

L 

Test section whole length 

[ml 

0 

input power 

rwi 

q 

Surface heat flux 

[W/m 2 1 

r 

Test section radial position measured from section center line 

[ml 

s 

Step height 

[ml 

T 

Temperature 

[K1 

U 

Local air velocity 

[ m/s 1 

V 

Heater voltage 

[voltl 

X 

Axial distance after step 

[ml 


Dimensionless Group 


Name 

Description 


Nu Dh 

Local Nusselt number based on hydraulic diameter 

(h x .Dh/k) 

Nu max 

Maximum Nusselt numberbased on hydraulic diameter 

( h x , max .Dh/k) 

Nu x 

Local Nusselt number based on axial distance 

( h x .X/k) 

Nu LDh 

Average Nusselt numberbased on hydraulics diameter 

( h[ .Dh/k) 

Reph 

Reynolds number based on hydraulic diameter 

(p.U.Dh/ p) 

Re x 

Local Reynolds number based on axial position 

(p.U.X/p) 

Re L 

Average Reynolds number 

(p L .U.Dh/ p L ) 

U* 

Dimensionless axial velocity 

(U/U.) 

R* 

Dimensionless radial position 

(2r-d)/Dh 
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Greek 


Symbols 

Description 

Units 

P 

Air film density 

( Kg/ur ) 

M 

Air film dynamic viscosity 

( Kg/m . s ) 

o 

Stephan-Boltzmann 

(W/mfk 4 ) 

8 

Surface emissivity value 



Subscript 


Subscript 

Description 

1,2 

Outer surface of inner tube and inner surface of outer tube average 

BX 

Local air bulk 

BL 

Average air bulk 

Cond. 

Conductive 

Conv. 

Convection 

FX 

Local air film 

FL 

Average air film for 1 length 

i 

Annulus average inlet (upstream) 

L 

Average for whole length 

max 

Maximum local value 

r 

Radiation 

sx 

Local surface 

SL 

Average surface 

X 

local 
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